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Employing a non-perturbative gauge invariant definition of the Debye screening mass mo in the 
effective field theory approach to finite T QCD, we use 3d lattice simulations to determine the leading 
0{g^) and to estimate the next-to- leading 0{g'^) corrections to mn in the high temperature region. 
The 0{g'^) correction is large and modifies qualitatively the standard power-counting hierarchy 
picture of correlation lengths in high temperature QCD. 
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QCD matter, a spatially and temporally extended sys- 
tem of matter described by the laws of Quantum Chro- 
modynamics, goes at high temperatures into a quark- 
gluon plasma phase, in which color is no more confined 
and chiral symmetry is restored. An essential quantity, 
describing coherent static interactions in the plasma, is 
the inverse screening length of color electric fields, the 
Debye mass mo- The Debye mass enters in many essen- 
tial characteristics of static properties of the plasma. Its 
numerical value is important for phenomenological dis- 
cussions of formation of the quark-gluon plasma, for the 
analysis of J/^! and T suppression in heavy ion colli- 
sions, for the computation of parton equilibration rates, 
etc. (see, e.g. 0). 

The definition and computation of the Debye mass for 
abelian QED plasma is well understood [||. The electro- 
magnetic current is a gauge-invariant quantity, and 
the Debye mass can be extracted from the 2-point gauge 
invariant correlation function of jo in the plasma. There 
are no massless charged particles in QED, which allows 
an infrared-safe perturbative computation of the Debye 
mass in powers of the electromagnetic coupling e. This 
has been done to order The situation in QCD is 

much more complicated. First, the corresponding cur- 
rent in QCD, is not a gauge invariant quantity. Sec- 
ond, there are massless charged gluons which give rise to 
infrared divergences and prevent the perturbative deter- 
mination of the Debye mass beyond leading order. 

A non-perturbative gauge invariant definition of the 
Debye mass in vectorlike theories with zero chemical po- 
tential was suggested in [Q. According to it, m/j can 
be defined from the large distance exponential fall-off of 
correlators of gauge-invariant time-reflection odd opera- 
tors O, 

(O(r,f)O(r,0)) - C\xfeM~mD\x\), (1) 
where C and /3 are some constants. The simplest choice 



for the operator O is Fq:^Fi2, and other examples can be 
found in In principle, 4-dimensional lattice simula- 
tions of hot QCD would thus allow a measurement of the 
Debye mass at any temperature. 

The aim of this letter is a non-perturbative determi- 
nation of the high temperature limit of the Debye mass, 
at T > a few x Tc. We will see that the effective 3d ap- 
proach to high temperature gauge theories, developed in 
(for a review, see |^) allows a simple and transpar- 
ent gauge invariant definition of the Debye mass |^ , while 
3d lattice Monte Carlo simulations provide an economical 
way to determine its value. The corrections to the leading 
result we shall find are numerically large; thus many com- 
putations in the phenomenology of quark-gluon plasma 
in heavy ion collisions should be re-analysed. 

The theory we shall study is QCD with Nf massless 
quark flavours and with the gauge group SU(iV) with 
iV = 2, 3. At high temperatures and zero chemical po- 
tential the Debye mass can be expanded in a power series 
in the QCD coupling constant g — g{fi) (the scale fi will 
be specified later; the result for Nf — is shown explic- 
itly in eq. (§)): 
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The leading order (LO) perturbative result, 
(7V/3 -|- Nf/Q^/'^gT, has been known for a long time . 
The logarithmic part of the 0{g^) correction can be ex- 
tracted perturbatively [p^ , but cat and the higher order 
corrections are non-perturbative. We are going to evalu- 
ate numerically the coefficients cjv and cIn^^^.. 

Static Green's functions for bosonic fields of high tem- 
perature QCD at distances ^ we are interested 
in can be determined by constructing an effective 3d 
gauge theory, containing static magnetic gluons and the 
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zero component of the 4d gauge field, Aq More- 
over, a super-renormalizable 3d theory, defined by the 
Lagrangian 



LMAl ^g] = -F^^F,^^ + Tr [A, Aq] [A, ^o] 
+ mlTTAl + XA{TTAl)^, 



(3) 



gives the Green's functions to a relative accuracy 
©(p"*) which is sufficient for the accuracy of the ex- 
pansion in eq. (||) . The parameters of the effective theory 
are related to the parameters of 4d QCD {A-^, N, Nf) 
and the temperature as described in | pH . For brevity, we 
give here the explicit expressions only for Nf — 0: 
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Here g{^) is the QCD coupling in the MS scheme and 
all the effective theory couplings have been computed 
including both the leading and the next-to-leading order 
contributions. The couplings (||)-(||) are independent of 
the gauge chosen for the perturbative computation. The 
expansion parameter is g"^ /l&i:^ so that the result should 
be accurate down to T r:! a few x Aj^. 

The dynamics of the 3d effective theory is fully char- 
acterised by the two dimensionless ratios 



mo 
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and by the dimensionful coupling g^. The value of x is 
essentially fixed by T, 
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while y and x, corresponding to physical 4d finite T QCD 
for Nf — 0, are related by: 

y = yM = ^ + ^ + 0{x), N^2 (9) 

We are now ready to give a gauge-invariant definition 
of the Debye mass in the 3d language Q . Physically, we 
want a local operator which makes Aq gauge invariant in 
the 3d theory and contains Aq singly. We can single out 
this state by a symmetry consideration. Note that the 
effective Lagrangian (||) has a discrete symmetry Aq ^ 
—Aq. Then the Debye mass can be defined as the mass of 



the lightest 3d state which is odd under this symmetry. 
The operator of lowest dimension of this type is 



hi — eijkT^i AqF, 



jk 



(11) 



Here the field Aq has been made gauge invariant by dress- 
ing it with a cloud of magnetic gluons. 

At high T, one has 5^1 and, according to eqs. (||), 
7713 ^ ffl- This is the "heavy quark hmit" of the 
3d theory, in which the mass niD of the singlet state is 
dominated by the bare mass of the scalar "quark" 
Aq. For dimensional reasons, the exact mass mo can in 
this limit be expanded as 



rriD = 7773 
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where un and 6jv are constants, perhaps involving a log- 
arithm of m^/g^. The terms neglected are of higher or- 
der using the power counting in eqs. (^-(§|). Comparing 
eqs. and (p^), one sees that 
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(13) 



Here we used the fact that the scale dependence of the 
non-perturbative terms in eq. (^ is at least of order 
0{g^). Since the expansion ([l^) refers only to the 3d the- 
ory, the constants and &jv depend on A^ but clearly 
not on Nf. Thus cn is A^^ independent, while d^^Nf 
depends on Nf only through . 

In terms of our dimensionless variables (^, eq. ( p^ ) 
becomes 



^ N ^ ^ bN 
47r y/y 



(14) 



The mass 777d can now be measured by putting the 
effective 3d theory on the lattice |l^ and by measuring 
the exponential falloff of the correlator (/73(a;3)/73(0)) ~ 
exp(— TMul^sl), where h3{x3) is summed over the trans- 
verse {xi, X2) plane. The correlation function is measured 
both with zero and finite transverse momentum, and in 
order to enhance the overlap with the asymptotic state 
the measurements are performed with several levels of 
recursive blocking of the operators. We select the block- 
ing level and momentum sector which has the best signal 
for the asymptotic mass separately for each Monte Carlo 
run. Since the longest correlation length in each case is 
less than 1 /5 of the linear size of the lattice, we expect the 
finite volume effects to be negligible in comparison with 
the statistical errors. This was also explicitly checked by 
performing simulations with different volumes in isolated 
cases. 
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FIG. 1. The gauge invariant Debye mass for SU(3), as a 
function of x, or T/Aj^ through eq. (H). The scale of g on the 
j/-axis is fixed according to eq. (^. The dashed line marks the 
leading order m^^, and the continuous line the 2-parameter 
fit to eq. (|l^) with the parameters as in eq. (|l5|). 



The mass is defined in the whole y,x parameter 
space. To have results which are relevant for 4d physics, 
we perform the measurements along the 2-loop dimen- 
sional reduction lines yd-[{x), eqs. (|-|10|). To measure 
the coefficients in eq. (|l^) one should use the part of this 
curve corresponding to ^ ^ 1. The results for = 3 
are shown in Fig. |l|, in units of 4d gT {— yj 3y /N in 3d 
units). The Monte Carlo runs are performed with sev- 
eral lattice spacings a, parametrised by /3g = 2N/{g'^a). 
For SU(3) /3g varies by more than an order of magnitude 
(although not at the same value of x), as shown in Fig. ||; 
for SU(2), the measurements are done with (3g — 20 and 
32. The top scale of Fig. |l| shows the physical temper- 
ature T/Aj^ along ydr(a;) -line. Note that the highest 
temperatures are larger than 10^*^° x Aj^ ^ lO^"*^ x Tc- 

At small X (large j/), the fit to the function ( p^ is 
very good, as indicated by the continuous line in Fig. |^. 
In order to see in detail the sensitivity of the fit to the 
parameters, in Fig. ^ we replot the SU(3) data (restricted 
to X < 0.05) in terms of the quantity 5m/ g'^ = rriD/g'^ — 
■y/y ~ ^luy^ as a function of 1/y^. The intersection 
of the curve with the vertical axis gives the value of cat 
and the slope gives 6jv = <ijv,JV/ \/^/3 -I- A//6. One can 
see that the linear fit is rather good even down to small 
values of ^/y. The large non-zero value of the intercept 
is very robustly determined. The slope &Ar is small and 
has a relatively large error. Only the statistical error is 
given, but the value of Bn also depends on the range of 
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FIG. 2. The corrections of order 0{gl) and 0{gl/^) to 
the SU(3) screening mass mo in 3d units, corresponding to 
eq. (111). 



y included. 

The result of the fits are 

SU(2): CAT = 1.58 ±0.20 ^at = -0.03 ± 0.25 

SU(3): CAT = 2.46 ±0.15 ^at = -0.49 ± 0.15 . 



(15) 



The large number cat is related to non-perturbative 3d 
effects, while the smaller dN,Nf can be viewed as being 

For A = 2 we 
is close to zero. 



related to the choice of scale in ^ 
can in practice only verify that d2.Nf 
Note that writing — Ncn, one has cat — 0.79 ± 0.10 
{N = 2), 0.82 ±0.05 (N = 3). 
One can observe the following: 

• The leading term is dominant only at extremely 
large T. For SU(3), the leading term is larger 
than the 0{g'^) correction for g < 1/2.46 or for 
r/Aj^>exp(87r22.46Vll)>10i^ This implies 
that the leading term only dominates when QCD 
anyway merges into a unified theory. 

• The four terms in eq. (^^ fit the data over all the 
range T > lOOTc rather well, and there is no need 
for further corrections. 



In the range w 7c < T< lOOTc, to_d is rather 
constant and « 3.0 for SU(2) and « 3.3 
for SU(3). It should be noted, though, that 
in this regime m]^ > T so that the hierarchy 
m|)/(27rr)^ <C 1 required for an accurate descrip- 
tion of 4d physics through a 3d effective theory is 
getting weaker. 
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• The mass measured from the (^o^o)-correlator m 
the Landau gauge in 4d simulations for Nf = has 
also been observed to be clearly larger than the 
leading term ||l2| . 

• If the mass ^ gT of the Aq field is "large" , larger 
than the non-perturbative 0{g^)T correction, the 
Aq field can be further perturbatively integrated 
out and a simpler effective theory, containing only 

(and possible scalar fields) can be derived 
Our results imply that this can be accurately car- 
ried out for QCD only at extremely high temper- 
atures, T ^ Tc- In the electroweak case the accu- 
racy of the integration is sufficient even for T ^ Tc 
both since the leading term has a bigger coefhcient 
(m^° = ^/n/6gT) than the Nf = QCD consid- 
ered here and since in the relevant T regime the 
coupling constant g = g{m\Y) ~ 2/3 is smaller. 

• The usual parametric "power counting" picture of 
correlation lengths in high temperature QCD says 
that the longest scale, related to the magnetic sec- 
tor of the theory, is m]^ ~ (const x g^T)^^. A 
shorter scale, ~ {gT)~^, is associated with Debye 
screening. Our results show that this picture can 
be quantitatively correct only at extremely large 
temperatures. Indeed, purely magnetic effects, as 
measured by the 3d glueball (operator F^jF^j) mass 
(mc « 2g| for pure SU(2) |T|,|n|l) tend to be nu- 
merically large, so that itim ^ mr, in a very wide 
range of temperatures (This gauge invariant result 
is in contrast to the small magnetic gluon masses 
measured in Landau gauge jlj]). In this range the 
longest length scale corresponds to a scalar 0++ 3d 
"bound state" of two Aq quanta, associated with 
the operator AqAq (the power counting suggestion 
that this state is roughly twice as heavy as mu 
holds again only at extremely high T). 

Summarizing, we have carried out with lattice Monte 
Carlo techniques a gauge independent measurement of 
the Debye mass in finite T QCD. The measurement is 
based on first deriving with 2-loop perturbative compu- 
tations a 3d effective theory. The expansion parameter 
is Qs/tt, so that the result is accurate down to T close to 
Tc- The mass is obtained by measuring correlators of the 
gauge invariant local operator AqF"^^. in the 3d theory. 
The leading and next-to-leading corrections to m^) were 
determined and found to be large. In fact, for tempera- 
tures from Tc up to T '--^ lOOOAj^ the non-perturbative 
Debye screening mass is about a factor 3 larger than the 
lowest order estimate. 

It remains to be seen whether this modification of 
the standard picture of high-temperature gauge theo- 
ries has applications in the cosmological discussion of the 
quark-hadron phase transition or in the phenomenology 
of heavy ion collisions. 



The simulations were carried out with a Cray C94 and 
Cray T3E at the Finnish Center for Scientific Computing. 
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